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ABSTRACT 



Aims. As a well-established procedure for the vast majority of normal main-sequence stars, determination of atmospheric and stellar 
parameters turns to be a challenging process in case of magnetic chemically peculiar stars. Inhomogeneous distribution of chemical 
elements and strong magnetic fields make most of the standard photometric and spectroscopic calibrations inapplicable for this class 
of stars. In this work we make use of available observed energy distributions calibrated to absolute units, stellar parallaxes, high- 
resolution spectroscopic observations, and advanced stellar atmosphere models to derive parameters of three bright Ap stars: 33 Lib, 
7 Equ, and f3 CrB. 

Methods. Model atmospheres and fluxes were computed with the LLmodels code. Synth3 and Synthmag codes were used to compute 
profiles of individual spectral lines involved in abundance analysis. 

Results. For each of the stars, we construct a self-consistent atmospheric models assuming normal and depleted helium compositions 
and derive empirically stratification profiles of certain elements. The effective temperatures and surface gravities are found from 
the simultaneous fit to spectroscopic, photometric, and spectrophotometric observations calibrated to absolute units. We show that 
using advanced model atmospheres and accurate stellar parallaxes allows one to derive stellar radii with high accuracy, and these are 
consistent with those obtained from independent but more complicated interferometric observations. 
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1. Introduction 

Accurate knowledge of fundamental stellar (mass, radius, lumi- 
nosity, metallicity, age) and atmospheric (effective temperature, 
surface gravity, atmospheric abundances) parameters is essential 
for many fields of modern astrophysics. Stellar masses and abun- 
dances determine the whole chain of stellar evolution stages of a 
single star, and obtaining stellar luminosities and radii by means 
of observed data allows one to determine ages of stars, as well as 
to derive other parameters of interest, such as atmospheric chem- 
istry from spectroscopic observations, to characterize extrasolar 
planets from transit data, to carry out asteroseismic modeling, 
and much more. 

Determining stellar parameters is a delicate process, even in 
the case of normal stars, because it relies on theoretical evolu- 
tionary and atmospheric models. The former requires knowing 
metallicities and masses, the latter also needs effective temper- 
atures and gravities as input. These are usually obtained from 
spectroscopic and photometric data employing well-developed 
calibration schemes. For chemically peculiar (CP) stars, how- 
ever, things get much more complicated mainly because of the 
abnormal chemical comp osition of their atmospheres (see, e.g., 
Ryabchikova et al. 2004). This makes the spectra and energy 
distributions of CP stars look very different from their normal 
analogs, and therefore most of the photometric calibrations sim- 
ply fail to work for CP stars. Furthermore, the strong surface 



Send offprint requests to: D. Shulyak, 
e-mail: denis . shulyak@gmail . com 



magnetic fields detected in many CP stars additionally compli- 
cate the analysis of their spectra and determination of atmo- 
spheric parameters. 

Since all parameters are linked, the accuracy in derived pa- 
rameters, such as T e ff and log(g), must be as high as possible. 
This is why the full set of stellar and atmospheric parameters 
can ideally only be derived using an iterative procedure where 
each parameter is adjusted in such a way as to fit all available 
observed data. Then, if needed, the other parameters are recal- 
culated accordingly. For instance, deriving T e g and log(g) from 
photometric data and employing relevant calibrations one can 
carry out a detailed abundance analysis based on model atmo- 
spheres, but if derived abundances turn out to be different from 
those assumed in the construction of the model atmospheres 
then the latter must be recalculated with the appropriate chem- 
istry, and the whole process should be repeated again. Because 
of the lack of appropriate calibrations and the wide diversity of 
abundance patterns, this is the only robust s cheme for analyzing 
CP stars (see e.g iKochukhov etai]l2009l: IShulvak et al] 2009: 
iPandev et alj|2011l) . 

Recently, interferometry has been successfully applied to 
CP stars resulting in a first estimate of the radii of three stars: 
a Cir (HP 12889 8, iBruntt et al.1 l2008h . B CrB (HD 137909 , 
iBrunttet al.ll2010l) . and y Equ (HD 201601. lPerraut et al.ll201 lb . 
No doubt interferometry is a powerful technique with the strong 
advantage of providing model-independent estimates of stellar 
sizes, but its current world-wide facilities, unfortunately, are lim- 
ited only to bright stars. 
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In this work we apply a self-consistent atmospheric analy- 
sis to three well-known magnetic CP stars: (5 CrB (HD 137909), 
y Equ (HD 201601), and 33 Lib (HD 137949). Collecting avail- 
able photometric, spectrophotometric, and spectroscopic obser- 
vations we perform detailed abundance and stratification analy- 
sis of a number of chemical elements in the atmosphere of each 
star. The effective temperatures and gravities are found itera- 
tively by comparing theoretical fits to all available observables. 
Using fluxes calculated consistently with model atmospheres, 
available stellar parallaxes, and observed fluxes calibrated to ab- 
solute units, we derive stellar radii and compare them with inter- 
ferometric results. 



2. Observations 

Our investigation is based mainly on the low-resolution spectro- 
scopic observations calibrated to absolute units and obtained by 
Space Telescope Imaging Spectrograph (STIS£] mounted at the 
Hubble Space Telescope. STIS provides wide wavelength cov- 
erage (/1/11700 - 10200) that is excellent for accurate analysis 
of the energy distributions of early type stars where both UV 
and visual fluxes are highly important. Among our targets, only 
33 Lib was not observed with STIS. Nevertheless, medium -band 
spectro-photometric observations by lAdelman et al.l d 19891) are 
available for all objects, as are UV fluxes from the International 
Ultraviolet Explorer (IUEfl For 33 Lib and ft CrB, additional 
wide -band UV fluxes were taken from the TD1 space mis- 
sion dBoksenberg et al.lll973l;lThompson et"aT1ll978l) . For/? CrB 
we also made u s e of the spectro-photometric measurements o f 
lAlekseevaetakl dl996l) . iBregerl (1 19761) . iKomarov et all (119951) . 
and UV low-reso lution fluxes from the Astron space satellite 
(lBoiarchukll986l) . 

We note that CP stars are known to demonstrate photo- 
metric variability, most of which is caused by the inhomoge- 
neous surface distributi on of chemical elem ents. However, as 
can be seen from e.g. [Shulyak et al] (l2010al) and Krticka et al. 
(2012), the amplitude of this variability is rather small to influ- 
ence the derived effective temperatures much. For instance, us- 
ing Adelman's spectrophotometric catalog, we find that the rel- 
ative characteristic correction to the effective temperatures due 
to flux variability is only a few ten K's for all three target stars. 
Because of this, and also because the datasets listed above are 
not homogeneous in the sense of observing times and number of 
observations performed, we used mean fluxes when possible. For 
instance, only phase-averaged fluxes were used from Adelman's 
catalog and the IUE data archive. 

Chemical abundance and stratification analysis was based 
on high-resolution, high sign al-to-noise-ra t io sp ectra obtained 
with the UVES instrument dDekker et al.l l2000h at the ESO 
VLT. Observations of 33 Lib and ft CrB were carried out 
in the context of progra m 68.D-0254, and published by, e.g., 
Rvabchikova et al. (2008), while spectra of y Equ were extracted 
from the ESO archive (program 76.D-0169). Spectra of all three 
stars were obtained using both dichroic modes, thus covering 
wavelength intervals from /13030 to /U0400 with a few gaps. 
Spectral resolution is about 80 000, and the signal-to-noise ra- 
tio at ,15000 is about 300- 400. A detailed de s criptio n of the 
data reduction is gi ven by Rvab chikova et al.l d2008l) (33 Lib 
and ft CrB) and by ICowlev et al.l (120071) (yEqu). In addition, 
for 33 Lib and y Equ we used the UVES sheer spectra cover- 
ing the /1/14960 - 6990 region and providing a resolving power 
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of A/AA ss 115 000. Time-resolved observation of these stars 
were obtained in the context of ESQ programs 72.D-0138 and 
077.D-0150 (e.g. iKurtz et al.l 120061) . We used average spectra 
constructed from the total of 111 and 178 short exposures for 
33 Lib and y Equ, respec tively. Reduction of these data followed 
procedures described by Rvabchiko va et al.l d2007l) . 

In the course of our analysis we found that for the j3 CrB 
and y Equ the H/3 line of the UVES spectra always require 
systematically different temperatures compared to the one de- 
rived from the spectral energy distribution (SED) and other hy- 
drogen lines. We attribute this to a problem of continuum nor- 
malization and made use of additional observations from differ- 
ent instruments. In particular, for/3 CrB we used spectroscopic 
data obtained at 2.6-m tele scope mounted at CrAO (Cri mean 
Astrophysical Observatory) dSavanov & Kochukhovi ri998). For 
y Equ we used medium resolution spectra (R =s 37 000) obtained 
with the echelle spectrograph GIRAFFE, mounted at the 1.9- 
m Radcliffe telescope of SAAO (South African Astronomical 
Observatory). The data reduction of the GIRAFFE spectra was 
based on Vienn a automatic pipeline fo r echelle spectra process- 
ing described in lTsvmbal et al.l d2003l) . 



3. Basic methods 

3.1. Model atmospheres 

Our approach is based on the models and synthetic fluxes 
computed with the L Lmqdels stellar model atmosphere code 
(Shul yak et al.l I2004I) . The code can treat individual homoge- 
neous and depth-dependent abundances, which is a key feature 
in our modeling. In addition, the effect of the magnetic field 
is fully accounted via detailed computations of the anomalous 
Zeeman splitting and pol arized radiative transfer as described in 
iKhan & Shulya k (2006a). In all calculations we adopted mag- 
netic field vect or to lay perpendicular t o the normal of the stel- 
lar surface (see IKhan & Shulyak 2006b, for more details). This 
makes LLmodels a powerful tool for the modeling of the ob- 
served characteristics of magnetic CP stars. The model fluxes 
are computed on a fine frequency grid with the highest point 
density at the wavelengths occupying the region of maximum 
stellar flux. A total of about 600 000 frequency points are usu- 
ally used, which provides the necessary resolution for accu- 
rate computation of photometric indicies of broad- and narrow- 
band photometric systems. Lines opacity coefficient is computed 
based on atomic parameters extracted fro m the VALD database 
dPiskunov et alj|1995t iKupka et aklll999l) . which presently in- 
cludes more than 6.6x 10 7 atomic lines. Most of them come from 
the latest theoretical calculations performed by R. Kuruc2Q. 

Because of their large overabundance in atmospheres of Ap 
stars, rare-earth elements (REE) play an important and some- 
times ev en dominating role in th e appe arance of spectroscopic 
features (Mashonkina et al. 2009, 2005) and can also affect the 
SED (Shu lvak et al.ll2010bl) . Here we used the same lin e lists 
of REE e lements as presented in Shulvak et al.l d2010bl) . The 
studies bv lMashonkina et all (I20091 12005b demonstrate that the 
line formation of Pr and Nd can strongly deviate from the lo- 
cal thermodynamic equilibrium (LTE) with the doubly ionized 
lines of these elements are unusually strong due to the combined 
effects of stratification of these elements and departures from 
LTE (known among spectroscopists as the "REE anomaly"). 
However, a precise NLTE analysis is beyond the scope of this 
paper and currently could not be directly coupled to a model at- 



3 http://kurucz.harvard.edu 



2 



D. Shulyak et al.: Fundamental parameters of bright Ap stars 



mosph ere calculation. There fore, we followed the approach out- 
lined in lShulvaket all (120091) where the authors used a simplified 
treatment of the REE NLTE opacity by reducing the oscillator 
strengths for the singly ionized REE lines by the corresponding 
LTE abundance difference, while using the abundances derived 
from second ions as input for model atmosphere calculations (or 
vice versa). Applying this line scaling procedure allowed us to 
mimic the line strengths that correspond to the NLTE ionization 
equilibrium abundance for each REE. 

3.2. Abundance and stratification analysis 

Initial homogeneous ab undances for mode l cal culations were 
taken from th e papers bvlRvabchikova et alJ d2004l) (33 Lib and 
(3 CrB) and by Rvabchik ova et al.l d 1997b for y Equ. These abun- 
dances were slightly corrected with synthetic spectral synthe- 
sis taking into account rather large magnetic fields of the pro- 
gram stars. Spectrum synthesis was perform ed by the magnetic 
synthesis code SYNTMAG (see lKochukhov 2007) using extrac- 
tion of the atomic parameters and Lande-factors from the VALD 
database. To make abundance analysis of magnetic stars more 
efficient, we used the equivalent width method with the newly 
developed code WidSyn. This program provides an interface to 
the SYNTMAG code and allows abundances to be derived us- 
ing theoretical equivalent widths obtained from the full polarized 
radiative transfer spectrum synthesis. This improved equivalent 
width method considerably accelerates abundance analysis with- 
out sacrificing accuracy. 

The stratification of chemical elements was modeled by 
employing the DDAFit IDL-based automatic procedure that 
finds ch emical abundance gradients from the observed spec- 
tra (see iKochu khov 2007). In this routine, the vertical abun- 
dance distribution of an element is described with the four pa- 
rameters: chemical abundance in the upper atmosphere, abun- 
dance in deep layers, the vertical position of abundance step, 
and the width of the transition region where chemical abundance 
changes between the two values. All four parameters can be 
modified simultaneously with the least-squares fitting procedure 
and based on observations of an unlimited number of spectral re- 
gions. This procedure was successfull y employed by a number 
of element stratification studies (e.g., [ Rvabchikova et al. 2005; 
IKochukhov et alJ2009HShulvak et alJ2009l:lPandev et alJ201 ll) . 

For Pr and Nd in y E qu, we performed an NLTE strat ification 
analysis as described in Mashonki na et al.l (120051 120091) using a 
trial-and-error method and the observed equivalent widths of the 
lines of the first and second ions. In NLTE calculations magnetic 
field effects were roughly approximated by pseudomicroturbu- 
lent velocity individual for each line. Such an analysis is appli- 
cable only for stars wi th small or moderate magnetic fields, e.g. 
for HP 24712. y Equ dMashonkina et al.ll2005l l2009h . 

3.3. Determination of fundamental parameters 

We follow the basic steps of the iterative proce dure of atmo- 
spheri c parameters determination outlined in, e.g JShuivak et al.l 
(2009). Briefly, this procedure consists of repeated steps of strat- 
ification and abundance analysis that provide an input for the 
calculations of model atmosphere until atmospheric parameters 
(r e fr, log(g), etc.) converge and theoretical observables (energy 
distribution, profiles of hydrogen lines, etc.) fit observations. 
Thus, the derived abundance pattern is consistent with the phys- 
ical parameters of the final stellar atmosphere. 



We used low-resolution spectroscopic and/or spectro- 
photometric data calibrated to absolute units and available stellar 
parallaxes to derive radii of investigated stars. Fitting the full en- 
ergy spectrum makes it possible to find reliable values for the 
atmospheric log(g) by matching the amplitudes of Balmer jump, 
which is very difficult or even impossible to do when only spec- 
troscopic indicators or photometric indices are used. The fit to 
observed energy distributions is carried out on each iteration of 
abundance analysis using a grid of model atmospheres computed 
for a set of [T e ff,log(g)] values. The radii of the star is found by 
minimizing the deviations between observed and predicted ab- 
solutely calibrated fluxes. 

However, one important uncertainty remains in this analysis. 
As follows from the computations of atomic diffu sion, helium al- 
ways sinks in subphotosheric layers of A-F stars (Michaud et al. 
1979). The same He depletion is als o found by the diffusio n 
models of the optically thin layers (iLeblanc & Moninl 12004). 
suggesting that the atmospheres of CP stars of spectral types 
A-F may well be He-deficient. No suitable diagnostic He lines 
are available for cool Ap stars to verify this conclusion obser- 
vationally. The question of whether helium is indeed removed 
from the atmosphere is particularly important because this may 
noticeably change the molecular weight and T-P stratification 
of He-weak models in deep layers where temperatures become 
high enough for He to contribute to the total opacity coeffi- 
cient. Having no clues to the true He stratification, we thus con- 
sider an additional case of He-weak atmosphere with adopted 
He abundance of [He/H]= -4 (decreasing He abundance below 
this value has only a marginal or no effect). 

4. Results 

In this section we present the results of the parameter determina- 
tions for individual Ap stars, which are combined in TableQ] For 
completeness, this table also lists our previous results derived for 
a few additional CP stars using the same analysis methodology. 
We also find that, with a few exceptions, the application of real- 
istic model atmospheres did not result in a significant change of 
the mean element abundance estimates with respect to the stud- 
ies cited above. Therefore, in the following we focus on the anal- 
ysis of the vertical abundance profiles. 

4.1./3CrB(HD 137909) 

Among the three targets, ft CrB is probably the most studied 
one. A number of independent spectrophotometric observations 
of energy distributions make it possible to test predictions of 
atmospheric models, as well as estimate uncertainties of de- 
rived fundamental parameters in detail. Furthermore, ft CrB is 
kn own to be a bina ry syste m with orbital el ements presented 
bv iTokovininl d!984l) and iNorth et al.l (I1998I) . More recently, 
Bruntt et alJ d2010l) have carried out interferometric observations 
of the system, which resulted in the first model-independent de- 
termination of the radii of primary (A) and secondary (B) com- 
ponents: R(A) = 2.63+0.09^ o , R(B) = 1 .56 + O.O7/? . However, 
authors still used scaled -solar models comp uted with ATLAS9 
model atmosphere code ( Kuruczl 19921 11993) to predict bolomet- 
ric flux of the system and thus to derive effective temperatures 
for both components. This was necessary because all modern 
flux data always contain the combined light of the two com- 
ponents. This is also true for the STIS observation, which was 
done with the 0.2" slit width on 25 Aug 2003 when the visible 
separation between system components was r < 0.1". Effective 
temperatures were found to be r e ft(A) = 7980 + 180/T and 
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T e ff(B) = 6750 + 230 K, respectively. The secondary component 
thus provides an important « 18% contribution to the bolomet- 
ric flux of the system and thus has been taken into account in our 
SED fitting procedure. 

Available spectra of B CrB made it possible to analyze the 
stratification of five chemical elements: Mg, Si, Ca, Cr, and Fe, 
using from 8 to 23 neutral and singly-ionized lines per element. 
We did not perform NLTE Pr-Nd stratification analysis because 
of a rather strong m agnetic field and an absen ce of any signifi- 
cant Pr-Nd anomaly dRvabchikova et aO2004l) . which is a man- 
ifestation of the rare-earth stratification. Starting from the T e s = 
8000 K, log(g) = 4. 3 homogeneous abundanc e model, and abun- 
dances presented in Rvabchikova et alJ (120041) . the atmospheric 
parameters were iterated to match the observed SED and pro- 
files of hydrogen lines. All model atmosphe res were computed 
with (B s ) =5.4 kG surface magnetic field (Rvabch ikova et ail 
2008). The best fitted parameter of the primary fall in the range 
8000 < T eS (A) < 8050 K, 3.9 < log(g)(A) < 4.0 and the respec- 
tive fits to the SED and hydrogen lines are shown on Figs.[T]and 
[2] Generally, the T e ff of the star is constrained by the slope of the 
Paschen continuum and hydrogen lines, while the log(g) is best 
controlled by the amplitude of the B aimer jump. We find that the 
contribution of the B CrB-B is crucial for the Paschen continuum 
and IR region, and when ignoring it we were not able to obtain 
a consistent fit to SED and H-lines with any set of atmospheric 
parameters. 

From Fig. [T] it is clear that there is a systematic difference 
between space and ground-based observations. The STIS, IUE, 
TD1, and Astron data agree very well in the UV region below 
/12900. At the same time, STIS provides less flux in the visual 
and NIR than all ground-based photometric and spectrophoto- 
metric observations, which in turn agree fairly well with each 
other, with the exception of NIR (where some deviation can be 
seen, especially in the case of data from the Odessa catalog) and 
/U3800 - 5000 region. Most likely this is due to the different 
calibration methods applied to the observed data, and it is hard 
to give a preference to any of these datasets. Because STIS pro- 
vides most homogeneous sets of observations covering a wide 
wavelength range from UV to NIR, we used it as a reference in 
the model-fitting process. 

Two problems still remain. First, none of the models can rea- 
sonably fit the UV spectrum of the star below A25QQ. The sec- 
ondary component is too faint to influence the combined flux 
in this spectral region. Because all UV data agree well with 
each other, we are inclined to conclude that this discrepancy be- 
tween theoretical and observed fluxes is due to missing opacity 
in our models. The second problem arises from the fit to H-lines. 
Only Ha and HB are fit satisfactory. The computed profiles of 
Hy are always wider than the observed ones even for He-weak 
model with lower T e ff = 8050 K that demonstrate somewhat bet- 
ter agreement for these two particular lines and worse for Ha. 
By increasing of the secondary within the limits given by 
Brunttetal. (2010), one can obtain a slightly lower T e g for the 
primary and thus narrower profiles of Hy, respectively, but the fit 
to both NIR fluxes and particularly Ha profile then becomes no- 
ticeably worse. By integrating STIS fluxes and theoretical fluxes 
of the A = [8050, 4.0, 2.50 + 0.07], B = [6750, 4.2, 1 .56] model, 
we estimated a corresponding uncertainty in the effective tem- 
perature due to a difference between observed and predicted 
fluxes to be AT^s = 50 K. This means that the true of B CrB 
is probably slightly overestimated in our analysis. On the other 
hand, an uncertainty in He abundance results in a similar 50 K 
error, thus accounting for both, these effects could lower a tem- 
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perature of the star by 100 K compared to the temperature de- 
rived in this study. 

The final adopted parameters are (component^ 
[T eS (K),log(g)(dex),R(R @ )]): A = [8050,4.0,2.50 + 0.07] 
for the He-weak model and A = [8100,3.9,2.47 ± 0.07 for 
the He-normal (e.g., with solar helium compositi on) model. 
The p arameters of the secondary were taken from iBruntt et al.l 
(120101) : B = [6750,4.2, 1.561. We used the Hippar cos parallax 
of 7T = 28.60 + 0.69 mas from Ivan Leeuwenl (120071) . The errors 
of the derived radii results from the stellar parallax uncertainty. 

Using ground-based spectrophotometry led to a noticeable 
increase i n the stellar r adii. For instance, the fit to the data ob- 
tained by Hregi3([l976|) requires A = [8050,4.0,2.70 + 0.07] for 
He-weak and A = [8050, 3.8, 2.65 + 0.07] for He-norm models, 
respectively. The differences between STIS and ground-based 
observations correspond to the non-negligible difference in the 
derived stellar radii of AR « Q.2R e . This reflects certain prob- 
lems in the applied flux calibration procedures. But in spite of 
the vertical offset, the shape of SEDs remains practically the 
same, therefore the effective temperatures derived from the two 
datasets are very similar. 

Figure [3] illustrates stratification profiles derived from the 
final He-norm and He-weak models. A substantial change in 
the position of the abundance jump is found for Mg, for which 
the He-norm model results in its narrowing and deeper position 
compared to the He-weak model. Elements Ca and Fe demon- 
strate changes in the abundance amplitude in upper atmospheric 
layers, but we find that solution for Ca is quite uncertain in the 
highest layers with estimated errors on the order of +1 dex or 
higher. Stratifications of Si and Cr stay almost unchanged, with 
Si accumulating at much deeper photospheric layers than what 
was found for other elements. The changes in stratification pro- 
files observed in Fig. [3] are a combined effect of decreasing He 
content and T e ff by 50 K, i.e. not entirely due to He content. 
The effect of decreasing He abundance manifests itself in the in- 
creased log(g) by 0.1 dex needed to fit the SED of the star. Thus, 
accounting for a possible peculiar He content does have a well 
determined effect on the stellar parameter determination, and we 
suggest to compute models with strongly deficient He content 
(i.e. in the same way it is done in this research) when analyz- 
ing cool CP stars where true He abundance cannot be derived 
spectroscopically. 
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Fig. 1. Observed and predicted energy distributions of /3 CrB for two best-fitted models with solar and depleted helium contents. 
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Fig. 2. Observed and predicted profiles of hydrogen lines of /3 CrB. Full line - He norm model; dashed line - He-weak model. 



4.2. y Equ (HD 201601) 

Another well known roAp star with a rich spectrophotomet- 
ric dataset available is y Eq u. Similar to B CrB , recent inter- 
ferometric observations by iPerraut et alj d201 lb provide a ra- 
dius of the star R = 2.20 + 0.\2R Q . It is also a binary system 
(see iFabricius et al.l l2002l) with y Equ-B having much cooler 
temperature (T eS (B) = [4570,4833] K) and thus smaller ra- 
dius than /3 CrB-B, which then translates t o the contribution o f 
about 6% to the total flux of the system dPerraut et al.l I201 lh . 
The authors also estimated the effective temperature of the sys- 
tem T eS (A + B) - 7364 + 235 K, and thus T eff {A) = 7253 K. 
This is much cooler than the temperature of T e ff = 7700 K 
adopted for y Equ in the latest detailed spectroscopic study by 
iRvabchikova et alJ d2002l) . 



Due to extremely sharp spectral lines and weaker mag- 
netic field which minimize the blending effects, y Equ pro- 
vides a possibility to perform a stratification analysis for 18 
chemical elements: Na, Mg, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, 
Co, Ni, Sr, Y, Zr, Ba, Pr, Nd. This represents the most com- 
plete empirical study of element diffusion in an atmosphere 
of a single star to date. Starting fr om the initial model wit h 
T eB = 7700 K, log(g) = 4.2 from lRvabchikova"et"aT1 ( fl997l) . 
we iterated stratification analysis and model atmosphere calcu- 
lations until a minimal deviation between observed spectropho- 
tometric observations and model fluxes was obtained. The sur- 
face magnetic field introduced in model atmosphere compu- 
tations was (B s ) = 4 kG. In Sc, V, Mn, Co stratification 
calculations laboratory data for hyperfine structure (hfs) were 
taken into account. They were extracted from the following 
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papers: IVillemoes et al.1 (Il992h (Sc I D. lArmstrong et all (1201 ll) 
(VII) iBlackwell- Whitehead et all (120051) (Mnl). iHoltetalJ 
d!999h (MnII), lPickering|(11996l) (CoD. lBergemann et al.ld2QloT) 
(Coll). Hyperfine structure for Pr n was also includ ed in the 
Pr stratification ana lysis (se e Mashon kina et al.ll2009l) with the 
/1/5-constants from Ginibre (1989). In the case of He-norm set- 
ting, we find that models with A = [7600,4.0,2.04 + 0.05], 
B = [4500 - 4750,4.5,0.65 - 0.70] provide the best fit to the 
observed SEP apply ing parallax n = 27.55 + 0.62 mas from 
Ivan Leeuwenl (I2007D . Assuming B = [5000,4.5,0.80] results in 
a higher temperature of the primary A = [7700, 4.0, 1 .95 + 0.05] . 
To fit the Balmer lines, however, a lower effective temperatures 
of about T e ff = 7300 K is needed. Therefore, to maintain a satis- 
factory fit to both Balmer lines and the SED, a slightly cooler 
A = [7550,4.0,2.07 ± 0.05] was chosen. We note that for- 
mally the model with log(g) = 3.8 provides a better fit to SED 
(X 2 = 1.19) than the model with log(g) = 4.0 Or" 2 = 1.40) with 
a difference visible at the Balmer jump region (top right plot of 
Fig. 3J. But the log(g) = 3.8 then results in too low a mass of 
the star, M - 0.98M o , which is unreasonable for a main se- 
quence object with T e ff = 7550 K. No distinction could be made 
between these two values of log(g) from hydrogen line profiles 
because they are insensitive to the atmospheric gravity in this 
temperature range. 

From the SED comparison presented in Fig. [4] it is seen that 
good agreement between observations and model predictions is 
obtained for all spectral intervals represented from the UV to IR. 
Ground-based spectrophotometric data are also found to be con- 
sistent with the STIS observations, except observations obtained 
by the Odessa research group which are systematically lower. 
The derived effec tive temperatures from the SED are reasonably 
close to those of IRvabchiko va et al.l ([1997). Thi s is because in 
the lat ter work the authors fitted observations of lAdelman et al.l 
( 1989), which required T e g = 7700 K to match the slope of the 
Paschen continuum. 

Figure |5]compares the observed and predicted profiles of hy- 
drogen lines. The normalization of UVES data at the HB line was 
found to suffer considerably from the data processing inaccura- 
cies. Therefore we used observed profile taken with GIRAFFE 
instrument at SAAO. The normalization of the STIS data in 
the same region is also not trivial but easier to perform due 
to its lower resolution and a wide wavelength coverage. On 
the other hand, Ha line of STIS data is still satisfactorly fit- 
ted with r e (j = 7550 K models, though the UVES profile lies 
systematically higher, and thus lower T e ff is needed to fit it. 
The underabundance of He formally resulted in an increase in 
log(g) from 3.8 to 4.0 so that the He-weak model parameters are 
A = [7550,4.0,2.06 + 0.05]. 

As follows from the results presented above, the most no- 
ticeably affected parameter is log(g), which underwent changes 
from 4.2 to 3 .8 -4.0 depending upon the helium abundance used. 
A val ue of log(g) = 4.2 was used in Rvab chikova et al.l (1 19971 
|2002|) to fit atomic lines using atmospheric models with a scaled 
solar abundances, while in our analysis we used model atmo- 
spheres and fluxes computed with realistic abundances of the 
star. Finally, the use of STIS data allowed us to achieve a more 
detailed fit to the region of the Balmer jump, whose amplitude is 
very sensitive to the value of log(g). 

The vertical distribution profiles of all 1 8 chemical elements 
derived with r e jf = 7550 K, log(g) = 3.8, He-norm model are 
presented in Fig. [6] All but REE elements (i.e. Pr and Nd) are 
found to accumulate at layers just above the photosphere. For 
elements Ti and Ba, stratification is not well constrained and 
probably is not present in the atmosphere of this star at all. A 
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Fig. 6. Depth-dependent distribution profiles of chemical ele- 
ments in the atmosphere of y Equ. 



cloud of Pr and Nd is observed at high atmospheres and the abun- 
dance gradients of these elements, together with Fe and Si, are 
the largest among the studied elements. 

4.3. 33 Lib (HD 137949) 

In the present investigation 33 Lib is the last target analyzed. 
It is a bright roAp star with many references available regard- 
less its spectroscopic and photometric properties. Unfortunately, 
no high-resolution spectroscopic energy distribution is available 
similar to what is provided by STIS, therefore a fit was done to 
the gr ound-based spectrophotometric data from Adelman et al.l 
d!989h . Our fitting began w i th T eS = 755 K, log(g) = 4.3 model 
used bylyabchikovaetaD (|203[2008). Stratification of Si, Fe, 
Ca, and Cr was derived and accounted for in the model atmo- 
sphere calculations. Average abundances of the REE are higher 
in 33 Lib than in y Equ by more than an o rder of magnitude. 
The s tronger magnetic field of (B s ) = 5 kG (IRvabchikova et all 
2008), together with high REE abundances, results in extremely 
strong blending that prevents the proper choice of spectral lines 
for any stratification study of other elements. It is also the rea- 
son we could not perform NLTE stratification analysis of Pr and 
Nd based on equivalent widths and pseudomicroturbulence ap- 
proximation of magnetic effects. It could only be done by prop- 
erly including the magnetic field in NLTE line formation code. 
Pr-Nd stratification is expected to be rather strong in 33 Lib 
because of the highest Pr-Nd anomaly observed in this star 
(Rvabchi kova et al.l2004l) : however, at present step we limit our- 
selves by taking the average observed REE anomalies into ac- 
count. 

The best-fit parameters are [7400,4.0,2.13 + 0.13] for both 
He-norm and He-weak models. As shown in Fig. [7] also 
[7300, 3.8, 2.19 + 0.13] and [7500, 4.0, 2.06 + 0.13] models pro- 
vide a reasonable fit, but with more pronounced deviations be- 
tween the observed and predicted SEDs at short and long wave- 
lengths. 

Similar to other cases described above, for 33 Lib we also 
found it difficult to fit the Balmer Ha line with the model param- 
eters required to fit the SED (see Fig . [8j . A higher T e ff is required 
to fit the observed profile of this line. On the other hand, HB and 
Hy lines are a satisfactory fit with the T e ff = 7400 K model. 

The depth-dependent distribution of elements derived us- 
ing the [7400, 4.0] He-norm and He- weak models is shown on 
Fig. [9] As expected, in both cases the stratification profiles are 
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Fig. 4. Observed and predicted energy distributions of y Equ for two best-fitted models with solar and depleted helium contents. 
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Fig. 5. Observed and predicted profiles of hydrogen lines of y Equ. Dash-dotted lines correspond to He-weak model. 



very similar, with large abundance jumps at deep atmospheric 
layers. 

5. Discussion and conclusions 

In an attempt to construct self-consistent models of Ap-star at- 
mospheres, we performed a simultaneous model fit to a num- 
ber of observed datasets, such as photometric wide-band fluxes, 
low- and middle-resolution spectrophotometric observations, 
and high-resolution spectra. Ideally, all of them must be fitted 
with a single model atmosphere and abundance pattern. The re- 
sults of the present investigation suggest that this is indeed pos- 
sible to accomplish by applying a realistic chemical composition 
and appropriate model atmospheres that account for the peculiar 
nature of CP stars. 



However, a few problems still exist. The most disturbing one 
is a systematic discrepancy between effective temperatures de- 
rived separately from SED and H-lines. In the case of /? CrB, 
only Ha and H/3 agree with SED, but Hy requires lower T e ff . For 
33 Lib the situation is the opposite with H/3 and Hy in agreement 
with the T e ff derived from SED, but Ho- clearly requires higher 
T e ff . A good match to SED, Ha, and H/3 was obtained for y Equ, 
though no high-resolution data was available for the Hy line re- 
gion. Since at least two out of three hydrogen lines always agree 
with the atmospheric parameters derived from fitting the SED, 
it is likely that we face the usual problem of continuum normal- 
ization of the hydrogen lines in cross-dispersed echelle obser- 
vations. Thus, one should always be cautious when determining 
atmospheric parameters solely based on the fit to hydrogen lines. 
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Fig. 7. Observed and predicted energy distributions of 33 Lib for two best-fitted models with solar and depleted helium contents. 




The abundance and stratification analysis was done using 
carefully selected spectral lines with accurately known transi- 
tion parameters, and yet there seems to be a missing absorption 
in the models of (3 CrB in UV region blueward A25Q0 (see top 
left panel of Fig.Q]). A secondary component is by far too faint 
and cool to influence the flux in that spectral region. 

For /3 CrB we find a noticeable deviation between the space- 
and ground-based observations of the SEDs in the visual and 
NIR regions. A natural reason could be a discrepancy between 
different flux calibration schemes, which by default are more 
challenging in the case of ground-based observations. This devi- 
ation is seen as a vertical offset in the observed fluxes. It does not 
affect the determination of (Fig. [I), but does affect the radius 
estimate by about 0.2R o . This is quite a noticeable uncertainty 



that highlights once again the need to have accurately calibrated 
fluxes for stellar parameter determinations. 

All previous spectroscopic studies resulted in the log(g) val- 
ues that were substantially over-estimated, as shown by direct 
fits presented on Figs.Q~||4] andQ Our values come mainly from 
the detailed fit to the amplitude and slope of the Balmer jump 
based on realistic model atmospheres. Previously, only H-lines, 
photometric co l ors and/or low-resolution spectrophotometry by 
Adel man et alJ (1 1 989t) were used in combination with scaled- 
solar metallicity models that did not account for individual abun- 
dances, magnetic fields, and improved REE opacity. 

We find that possible He depletion of Ap-star atmospheres 
has minor effect on determining element distributions and log(g), 
at least for the hottest star considered here, /3 CrB(see Fig. 0. 
Formally, an underabundance of He results in choosing a slightly 



8 



D. Shulyak et al.: Fundamental parameters of bright Ap stars 



-2 




-5 -4 -3 -2 -1 

logr R 





-2 
-4 
-6 


1 1 1 1 

■ + s: 

■ % Co 

. O Cr 

■ A Fe 


1 1 1 1 1 1 1 1 

t7400g4.0 -»t» ■ 

t7500g4.0 

jr. ■ ■jytTt 1 1 

/ M ff ~ 




-8 




1 




ff * L 1 ' 

ft ,' / 1 ; 


log ( 


-10 




It /J / - 




-12 


^ / V 










-14 










-5 -4 -3 -2 -1 

logr R 



Fig. 9. Depth-dependent distribution profiles of chemical elements in the atmosphere of 33 Lib for He-norm and He-weak = 
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Table 1. Fudamental parameters of Ap stars. 



Name 


Zdf,K 


l°g(£) 


He-norm 
R, R Q L, L q 


M, M B 


7W> K 


log(j) 


He-weak 
R, R L, L 


M, M Q 


^inler' ^® 


(B s ), kG 


HD 137909 


8100 


3.9 


2.47 ± 0.07 


23.69 ± 1.93 


1.77 ±0.51 


8050 


4.0 


2.50 ± 0.07 


23.67 ± 1.91 


2.28 ± 0.65 


2.63 ± 0.09 1 


5.4 


HD 201601 


7550 


4.0 


2.07 ± 0.05 


12.56 ±0.94 


1.56 ±0.44 


7550 


4.0 


2.06 ± 0.05 


12.44 ±0.93 


1.55 ±0.43 


2.20 ±0.12- 


4.0 


HD 137949 


7400 


4.0 


2.13 ±0.13 


12.27 ± 1.83 


1.66 ±0.58 


7400 


4.0 


2.13 ±0.13 


12.27 ± 1.83 


1.66 ±0.58 


5.0 


HD 24712-' 


7250 


4.1 


1.77 ±0.04 


7.81 ±0.57 


1.44 ±0.40 




3.1 


HD 101065" 


6400 


4.2 


1.98 ±0.03 


5.92 ± 0.37 


2.27 ± 0.59 




2.3 


HD 103498 5 


9300 


3.5 


4.56 ± 0.77 


140.28 ±50.39 


2.40 ± 1.36 
















9500 


3.6 


4.39 ± 0.75 


141.57 ±51.35 


2.80 ± 1.60 
















HD 128898" 


7500 


4.1 


1.94 ±0.004 


10.74 ±0.33 


1.73 ±0.41 




1.967 ±0.066' 


~ 2.0 



The error bars of M and R are average errors computed taking parallax uncertainly and additionally assuming errors of AT^g 
Alog(g) = 0.1 dex into account. Larger errors found for HD 103498 result from a large parallax uncertainity. 
a Radius derived by me ans of interferometr y. 



50 K and 



Bruntt et al. 



Bruntt et al. 



2010); 2 lPerrautetal] ( l20TTI) : 3 IShulvak et td1j2009h : 4 IShulvak et al] ( l2010bl) ; 5 lPandev et al]d20Tlh : 6 KochuM ioTetal] (12009); 



2008) 



higher log(g) by about 0.1 dex and a cooler T e g by as 50 K. This 
change is small and well within the error bars of our fitting pro- 
cedure. The choice of temperature has a somewhat stronger or 
at least a comparable impact on the derived stratification profiles 
of chemical elements, as can be seen from the example of 33 Lib 
illustrated in left- and right-hand panels of Fig. [9] In particular, 
a change of Ar e ff = 100 K influences both positions and am- 
plitudes of the abundance jumps. It is important to realize, how- 
ever, that even if we are not absolutely precise in deriving r er r 
and know nothing about the real He content of the Ap-star atmo- 
spheres, the inferred element stratifications are still informative 
and stable. Additional uncertainties could result from particular 
limitations in our fitting procedure (e.g., step-like approximation 
of stratification profiles, a small number of spectral lines used, 
etc.), and thus the presented analysis can only provide a gen- 
eral picture of the element distributions in the atmospheres of 
investigated stars. Still, this is enough to reveal an accumulation 
of a particular element in the lower or upper atmospheric lay- 
ers, determine the amplitudes of abundance jumps, and estimate 
uncertainties in the derived distributions due to assumed model 
parameters - questions we were mainly interested in. 

Interferometry has recently brought an entirely new possibil- 
ity of measuring the radii of main-sequence stars on a model- 
independent and regular basis, though its modern instrumen- 
tal capabilities are limited to bright nearby stars. Our model- 
based comparison of the observed and predicted SEDs led to 



good agreement (within error bars) with interferometric esti- 
mates as summarized in Table \T\ The table also includes cal- 
culated masses and luminosities from derived effective tempera- 
tures, radii, and gravities. The corresponding error bars were cal- 
culated using individual uncertainties of stellar parallaxes along 
with typical uncertainties of Ar er y = 50 K and A log(g) = 0.1 dex 
associated with the fit of SEDs and hydrogen line profiles. Large 
error bars found for HD 103498 are entirely due to a large par- 
allax u ncertainty (n - 3.37 + 0.56 mas given by Ivan Leeuwenl 
rt2007k 

Figure[10]shows the positions of stars from TableQ]in the H- 
R diagram along with theoretical evolutionary tracks. This com- 
parison allows deriving stellar masses and ages under a reason- 
able assumption about the bulk stellar metallicity. Our spectro- 
scopically derived masses generally show good agreement with 
the theoretical values predicted by evolutionary models within 
estimated error bars. The only exception is HD 10 1065 for which 
evoluti onary models predict M — 1 .48 M , while Shulva k~et al.l 
(2010b) estimates 2.27 ± 0.59 M . In this particular case, how- 
ever, the derived log(g) was very sensitive to the assumed line 
opacity of REE elements that were found to affect model fluxes 
in the B aimer jump region in a dramatic way. Because of NLTE 
effects and the possible incompleteness of REE line lists, the 
published value of log(g) = 4.2 derived from the fit to the 
observed SED could easily be overestimated, and a somewhat 
lower log(g) = 4.0 would give a consistent mass estimate. 
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Fig. 10. Predicted location of Ap stars in the H-R diagram for the 
luminosity and temperature estimates from Table Q] Theoretical 
evolutionary tracks were adopted from iGirardi et al.l d2000l) for 
the bulk metallically of Z — 0.018. The dotted lines show 
ZAMS, 50% fractional main sequence age line, and TAMS. 
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